Vacuum Fluctuations of Energy Density can lead to the observed Cosmological 

Constant 
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The energy density associated with Planck length is puv oc Lp^ while the energy density associated 
with the Hubble length is pir oc LJj where Lh ~ 1/H. The observed value of the dark energy density 
is quite different from either of these and is close to the geometric mean of the two: pvac — ^puvpir- 
It is argued that classical gravity is actually a probe of the vacuum fluctuations of energy density, 
rather than the energy density itself. While the globally defined ground state, being an eigenstate 
of Hamiltonian, will not have any fluctuations, the ground state energy in the finite region of space 
bounded by the cosmic horizon will exhibit fiuctuations Apva,c{Lp,LH)- When used as a source of 
gravity, this Ap should lead to a spacetime with a horizon size Lh- This bootstrapping condition 
leads naturally to an effective dark energy density Ap ex {LuvLh)~^ oc H^ /G which is precisely the 
observed value. The model requires, either (i) a stochastic fluctuations of vacuum energy which is 
correlated over about a Hubble time or (ii) a semi-anthropic interpretation. The implications are 
discussed. 

PACS numbers: 



The conventional discussion of the relation between 
cosmological constant and vacuum energy density is 
based on evaluating the zero point energy of quantum 
fields with an ultraviolet cutoff and using the result as 
a source of gravity. Any reasonable cutoff will lead to a 
vacuum energy density pvac which is unacceptably high 

Q 

This argument, however, is too simplistic since the zero 
point energy — obtained by summing over the {l/2)hujk 
— has no observable consequence in any other phenom- 
ena and can be subtracted out by redefining the Hamil- 
tonian. The observed non trivial features of the vacuum 
state of QED, for example, arise from the fluctuations (or 
modifications) of this vacuum energy rather than the vac- 
uum energy itself. This was, in fact, known fairly early 
in the history of cosmological constant problem and, in 
fact, is stressed by Zeldovich ^ who explicitly calculated 
one possible contribution to fluctuations after subtract- 
ing away the mean value. This suggests that we should 
consider the fluctuations in the vacuum energy density 
in addressing the cosmological constant problem. 

Similar viewpoint arises, more formally, when we study 
the question of detecting the energy density using grav- 
itational field as a probe. Recall that an Unruh-DeWitt 
detector with a local coupling Lj — M{t)(J)[x{t)] to the 
field 4> actually responds to (O|0(a;)0(?/)|O) rather than 
to the field itself Q- Similarly, one can use the gravi- 
tational field as a natural "detector" of energy momen- 
tum tensor Tab with the standard coupling L = nhatT"'''. 
Such a model was analysed in detail in ref. Q and it 
was shown that the gravitational field responds to the 
two point function {0\Tab{x)Tcd{y)\0) . In fact, it is es- 
sentially this fluctuations in the energy density which is 
computed in the inflationary models |5j as the seed source 
for gravitational field, as stressed in ref. [g. AH these 
suggest treating the energy fiuctuations as the physical 
quantity "detected" by gravity, when one needs to incor- 



porate quantum effects. 

This fact alone, however, would not have helped. The 
vacuum expectation values in question are all ultravio- 
let divergent and will — normally — scale as a suitable 
power of the UV-cutoff, Luv If the energy density be- 
haves as p ex L^^, the fluctuation Ap will also scale as 
L~^ since that is the dominant scale in the problem; so 
the cosmological constant will still come out too large. 
Another difficulty is that, at least formally, the ground 
state is an energy eigenstate and it will have no disper- 
sion in the energy. However, as we shall describe, there is 
a nice manner in which De-Sitter spacetime circumvents 
these difficulties and leads to the correct value for the 
cosmological constant. 

The key new ingredient arises from the fact that the 
properties of the vacuum state (or, for that matter, any 
quantum state in field theory) depends on the scale at 
which it is probed. It is not appropriate to ask ques- 
tions without specifying this scale. (This, in some sense, 
has been the key lesson from renormalisation group.) If 
the spacetime has a cosmological horizon which blocks 
information, the natural scale is provided by the size of 
the horizon, Lh, and we should use observables defined 
within the accessible region. The operator H{< Lh), 
corresponding to the total energy inside a region bounded 
by a cosmological horizon, will exhibit fluctuations AE 
since vacuum state is not an eigenstate of this opera- 
tor. The corresponding fluctuations in the energy den- 
sity, Ap oc (A£')/L|^ = /{LuvtLh) will now depend on 
both the ultraviolet cutoff L„„ as well as Lh- We now see 
an interesting possibility of boot strapping: When used 
as the source of gravity, this Ap should lead to a space- 
time with the horizon size Lh, which, in turn, requires 
us to compute Ap using Lh as the infrared cutoff scale. 
Note that it is the existence of a cosmological horizon 
at Lh which provides a clear justification for using this 
length scale in computing the energy fluctuations. This 



leads to a value for Ap different from L~^ since we now 
have two length scales in the problem. This bootstrap- 
ping will lead to a relation between Ap, L^v and Lh- We 
will show that this relation is Ap oc [LuvLh)^"^ oc H^/G 
if we take the ultraviolet cut off at Planck length. That 
is precisely what we need. 

Remarkably enough, this result is very easy to obtain. 
To obtain Apvac oc AE/L^ which scales as (LpLh)^^ 
we need to have {AEy oc Lp L^; that is, the square of 
the energy fluctuations should scale as the surface area 
of the bounding surface which is provided by the cos- 
mic horizon. A simple argument to show that such a 
scaling is likely to occur is the following. Let the total 
Hamiltonian H of the system be written a.s H = Hi+ H2 
where Hi is the Hamiltonian obtained by integrating the 
Hamiltonian density within a sphere of radius R{= Lh) 
and 7?2 denoting the Hamiltonian of the outside region. 
(We are dividing the real space rather than separating 
scales in the momentum space, as is done in standard 
RG analysis.) Let Ei and E2 be the expectation values 
of Hi and H2 in the ground state |0) which is taken to 
be an eigen state of H with H\Q) — E\Q). Consider now 
the energy dispersions in the regions 1 and 2. Writing 
(ifa - E2f as [{H - E) - {Hi - Ei)]^ and noticing that 
the expectation values of {H — E) and {H — E)'^ vanish 
in any energy eigenstate, we get the result 

{AE2f = mH2~E2f\0) = {0\iHi-Eif\0) = [AEif 

(1) 
That is, the dispersions in the energy in the regions 1 
and 2 are equal. Since the regions 1 and 2 only share 
the bounding surface, this relation suggests that either 
dispersion could be proportional to the area of the surface 
and will scale as R^. The dependence on Lp is then 
fixed by dimensional considerations and we get (AE)'^ ex 
Lp R . 

While the above argument is suggestive, it leaves sig- 
nificant scope for improvement — not in the least because 
one is dealing with formally divergent quantities. Fortu- 
nately, it is possible to do this calculation rigorously. In 
flat spacetime, a scalar fleld has the familiar mode ex- 
pansion 



has the scaling 
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in terms of the harmonic oscillator modes qk(i)- Express- 
ing the Hamiltonian in the region 1 



Hi = 



x\<Lf 



d3x-(<^2-f|(V</.)p) (3) 



in terms of q\^ and gk one can evaluate the dispersion 

{AEf ^ {0\{Hi - {Hi))'\0) (4) 

in a straightforward manner. To obtain a flnite result, 
one needs to use an ultra violet regulator which we take 
to be the Planck length. For Lh ^ Lp, the flnal result 
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where the constant ci depends on the manner in which 
ultra violet cutoff is imposed. Similar calculations have 
been done (with a completely different motivation, in the 
context of entanglement entropy) by several people and 
it is known that the area scaling found in Eq. Q , propor- 
tional to i|j, is a generic feature [3. For a simple expo- 
nential UV-cutoff, ci = (I/SOtt^) but we do not believe 
this can be computed reliably without knowing the full 
theory. The result should also scale with the total num- 
ber of degrees of freedom at sub-Planckian energies which 
is not reliably known. [It is amusing to note that, the 
standard Zeldovich-Harrison scale invariant spectrum has 
the mass fluctuation {SM/Mf oc k^P{k) oc fc'' oc R-^ 
leading to the same "area" scaling (SM)'^ oc R^.] 

We thus find that the fiuctuations in the energy density 
of the vacuum in a sphere of radius Lh is given by 



AE ,_2.-2 H^ 
Apvac = ^3- OC Lp L^ oc 



G 



(6) 



The numerical coefficient will depend on ci as well as the 
precise nature of infrared cutoff radius (like whether it is 
Lh or Ln/'in etc.). It would be pretentious to cook up 
the factors to obtain the observed value for dark energy 
density. 

The current result in which the fluctuations in the en- 
ergy density act as a source of gravity is similar in spirit 
to the idea presented earfier in ref.|3]- If a quantum 
system is in a stationary state, its phase will evolve as 
exTp{—iEt). But for a partial system which can undergo 
fluctuations in the energy, one cannot attribute such a 
phase and we will — instead — obtain an uncertainty re- 
lation of the kind AE oc (1/t) where r is the time scale 
of the fluctuations. Similarly, the semiclassical gravity 
will have a wave function with a phase exp(— ip^acV) in 
case of a constant vacuum energy density, where V is 
the four volume. In case of a finite region of space, the 
fiuctuations lead to the result Apvac oc {LhLp)~^. To 
be rigorous, we need to do this computation in the De- 
Sitter geometry rather than in fiat spacetime. There is no 
conceptual difficulty in carrying out this program which 
merely requires replacing the mode functions in Eq. ^ 
by the ones appropriate for De-Sitter spacetime. It is 
easy to see from dimensional considerations that the re- 
sults in Eq. (Q and Eq. 10 will continue to hold. The 
numerical factor changes and needs to be determined by 
numerical integration. 

We have concentrated on the energy density rather 
than on the full energy momentum tensor for the sake 
of simplicity. More formally, one can construct the fluc- 
tuations in Tab along the same lines thereby obtaining the 
equation of state for the dark energy. If the regulariza- 
tion of the divergent expressions is handled in a Lorentz 
invariant manner (like, for example, using dimensional 



regularization in curved spacetime) , then the Lorentz in- 
variance dictates that the final result should have the 

form ATab = Apvacffab- 

The current observations ( 9] ; for earlier indications of 
nonzero cosmological constant, see Jj]j ) do suggest that 



Pvac — ^/PuvPiR 



(7) 



if we take the UV scale as Planck scale and the IR scale 
as Hubble scale. This remarkable relation deserves an 
explanation which is provided in a simple manner by our 
analysis. In fact, if we start with this result and do a 
bit of "reverse-engineering" , we immediately obtain the 
area scaling of energy fluctuations. The latter has been 
pursued and obtained by several authors, in connection 
with the holographic ideas and entropy-area connection. 

The main criticism one could raise is regarding our as- 
sumption that we can ignore the mean energy density, 
regulated by the UV cut off to a large value and con- 
centrate on its fluctuation. In addition to the arguments 
given in the first paragraph of this paper, there is an- 
other aspect which need to be stressed: It is a fact of 
life that a fluctuation of magnitude Ap^ac — H^ /G will 
exist in the energy density inside a sphere of radius H^^ 
if Planck length is the UV cut off. One cannot get away 
from it. On the other hand, observations suggest that 
there is a p^ac of similar magnitude in the universe. It 
seems natural to identify the two, after subtracting out 
the mean value by hand. Our goal was more towards ex- 
plaining why there is a surviving cosmological constant 
which satisfies Eq.Q which — in our opinion — is the 
problem. Another technical objection could be regarding 
our computing the /S.E over a sphere of radius H~^ and 
dividing it by the volume to get Ap^jac- This corresponds 
to assuming the coherence scale of fluctuation to be Hub- 
ble radius, which is indeed the correct thing to do. Recall 
that it is the horizon at H^^ which motivates choosing 
this scale rather than some other smaller scale. 

How does one interpret this result ? There are two 
possible ways. The first is to claim that the currently 
observed value of cosmological constant has no special 
significance and the fluctuations will lead to a p^ac ~ 
H'^{t)/G at any epoch. This suggests a stochastic inter- 
pretation of cosmological constant. To obtain the evolu- 
tion of the universe, we now need to solve the Friedmann 
equations with p^ac being a stochastic variable with a 
probability distribution V{pvac) which can also be com- 
puted in the ground state. Since the energy fluctuations 
can be positive or negative, in general, the cosmological 
constant can have either sign. There are (at least!) two 
difficulties in this interpretation: 

(i) For this idea to be viable, we need the fluctua- 
tions to be correlated over a time scale of the order of 
H(t)~^. It is trivial to compute the two point func- 
tion C{t) = {0\Hi{ti)Hx{tx + t)\0) for the quantum fluc- 
tuations and one flnds that C{t) = C2{L'jjt~^); this 
quantity is finite and we do not need to introduce the 
UV cutoff. (If the cutoff is introduced the result is 



C{t) ex L|j(Lp -I- t^)~^). Being a power law correla- 
tion, this suggests fluctuations at all time scales with 
the scale governed by Lh- But when t « Lh, the ratio 
[C(iH)/C(0)]i/2 = {Lp/Luf < 1 ifCifi) is computed 
with a cut off.. This probably means that coherence over 
t w Lh is rare. But the argument is not conclusive since 
what is required is the coherence time, given that a fluctu- 
ation of order L^/Lp has occurred. If we assume that co- 
herence over t « Lh is possible, then the current positive 
value of cosmological constant is a result of a fluctuation 
lasting typically for a time scale of Hq and we would 
predict to see a reversal of sign of pvac at earlier epochs. 
(Except that, it is difficult to predict whether this should 
have occurred at z = 2 or at z = 27r'^, say, with any con- 
fidence until all the numerical factors are accounted for. 
On the other hand, if the SN data at a redshift of 2, 
say, does indicate a negative cosmological constant, ones 
faith in this stochastic fluctuations will increase!). With 
this interpretation, the model solves both the conven- 
tional "problems" of the cosmological constant. Its value 
is correctly explained as due to energy fluctuations and 
this relation holds at any epoch over a local time scale of 
H{t)~^ thereby answering "why now?" by "always!". 

(ii) The second difficulty is is not special to our model 
and arises whenever one tries to model quantum fluctu- 
ations as a classical stochastic process and treat it as a 
source of semiclassical gravity. In the standard cosmo- 
logical model, the consistency of ([]) component of Ein- 
stein's equation with the (J) component will demand that 
d{pa^) — —pd{a^) which in turn will require p to be in- 
dependent of time ii p = —p. This relation will require 
a reinterpretation when cosmological constant is treated 
as a stochastic variable. If the coherence time of the 
fluctuation is of the order of Hubble time, there will be 
constancy of effective pvac (and standard evolution of the 
universe) within one Hubble time. But when stochastic 
fluctuation of cosmological constant occurs one needs to 
interpret the energy conservation equation in some lim- 
iting form. 

Other issues that one could raise regarding the com- 
putation in this approach are relatively straight forward 
to address: one can work out the probability distribu- 
tion for fluctuations in Tab in a De-Sitter background 
and integrate the stochastic Friedman equation. One 
also need to study the effect of the cosmological constant 
fluctuations on the growth of structures. All of these are 
best addressed numerically since one need to deal with a 
stochastic process with flnite correlation. 

There is a completely different way of interpreting 
this result based on some imaginative ideas suggested by 
Bjorken |ll| recently. The key idea here is to parametrise 
the universes by the value of Lh which they have. It is 
a flxed, pure number for each universe in an ensemble 
of universes but all the other parameters of the physics 
are correlated with Lh- This is motivated by a series of 
arguments in ref. [Tl| and, in this approach, p^ac oc L~^ 
almost by definition; the hard work was in determining 
how other parameters scale with Lh- In the approach 



suggested here, a dynamical interpretation of the scahng 
Pyac oc L^ is given as due to vacuum fluctuations of 
other fields. We now reinterpret each member of of the 
ensemble of universes as having zero energy density for 
vacuum (as any decent vacuum should have) but the ef- 
fective pvac arises from the quantum fluctuations with the 
correct scaling. One can then invoke standard anthropic- 
like arguments (but with very significant differences as 
stressed in ref. jlj ) to choose a range for the size of our 
universe. This appears to be much more attractive way of 
interpreting the result than introducing stochastic, time 
dependent fluctuations. 



Finally, to be fair, this attempt should be judged in 
the backdrop of other suggested solutions almost all of 
which require introducing extra degrees of freedom in 
the form of scalar fields, modifying gravity or introducing 
higher dimensions etc. and fine tuning the potentials (for 
a non-comprehensive sample of references, see J12j.^ At a 
fundamental level such approaches are unlikely to provide 
the final solution. 

I thank K.Subramanian for several rounds of extensive 
discussions and R.Nityananda and A.D.Patel for useful 
comments on previous drafts. 
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